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Abstract — A multivariate interpolation formula (MVIF) over 
finite fields is presented by using the proposed Kronecker delta 
function. The MVIF can be applied to yield polynomial relations 
over the base field among homogeneous symmetric rational 
functions. Besides the property that all the coefficients are 
coming from the base field, there is also a significant one on the 
degrees of the obtained polynomial; namely, the degree of each 
term satisfies certain condition. Next, for any cyclic codes the 
unknown syndrome representation can also be provided by the 
proposed MVIF and also has the same properties. By applying the 
unknown syndrome representation and the Berlekamp-Massey 
algorithm, one-step decoding algorithms can be developed to 
determine the error locator polynomials for arbitrary cyclic 
codes. 

Index Terms — cyclic codes, homogeneous symmetric rational 
function, Kronecker delta function, multivariate interpolation 
formula, unknown syndrome representation. 



I. Introduction 



A. Earlier Work 



Cyclic codes, proposed by Prange in 1957 1 1|, are important 
and practical error correcting codes. They are widely used 
today, including the Bose-Chaudhuri-Hocquenghem (BCH) 
codes and Reed-Solomon (RS) codes. Many decoders have 
been developed since then, including the famous Berlekamp- 
Massey (BM) algorithm which was designed specifically for 
the BCH codes or RS codes at first. Actually, other cyclic 
codes can also be decoded by the well developed BM al- 
gorithm, provided that there are enough consecutive known 
syndromes; namely, 2i consecutive syndromes are needed to 
correct a corrupted word with at most t errors. Unfortunately, 
for any cyclic codes other than the BCH/RS codes, the number 
of consecutive known syndromes is less than 2t. To obtain the 
unknown syndromes, Feng and Tzeng [2] proposed a matrix 
method which yields expressions of unknown syndromes in 
terms of known syndromes. With the extra syndromes, some 
BCH codes can be decoded up to their true error correcting 
capacity rather than the smaller designed one. Later, He et al. 
O developed a modified version of Feng's method, and used 
it to determine the needed unknown syndrome and then to 
decode the binary quadratic residue (QR) code of length 47. 
By applying the Feng-He's matrix method accompanied by 
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the BM algorithm, algebraic decoders are presented, by the 
authors, for the binary QR codes of lengths 71, 79, 89, 97, 
103, and 113, respectively 0-Q. 

A strict requirement to apply the original Feng's or the 
modified He's algorithms is that the weight of error pattern 
must be given. This leads to step-by-step decoding algorithms 
and then the error locator polynomial may not be determined 
in one step. Moreover, when the code length increases, the 
desired matrices do not exist in high-weight error cases. For 
example, no single matrix can be used in the weight-5 error 
case of the QR code with length 47; actually, He et al. used 
two matrices rather than one. To decode each of the six QR 
codes of lengths 71 through 113, the authors Q-Q used at 
least two matrices. 

In order to develop one-step decoders, the Lagrange interpo- 
lation formula (LIF) was applied [4| to yield representations of 
the unknown syndromes in terms of known syndromes when 
the codes are generated by irreducible polynomials. The LIF 
is a well known technique. For any finite set of points in the 
xy-plane with all the x-coordinates distinct, the LIF can be 
used to provide a polynomial function whose graph passing 
through these given points. The LIF can be applied in finite 
fields. And, in the unknown syndrome representation case, the 
obtained polynomial has nice properties: All its coefficients are 
coming from the base field; moreover, the degree of each term 
has the same remainder when divided by the code length, as 
mentioned in the following theorem appeared in |4|. 

Theorem A. (LIF) Let Q = {e(/3) | 1 < wt{e{x)) < t} be 
the set of all correctable syndromes with t the error correcting 

capacity. For a = (3 h + h ft lv efl with < h < ■ • • < 

l v < n - 1, 1 < v < t, let aM = /3 rh H h P rU with 

1 < r < n — 1. Then the polynomial L r (x) defined below 



L r (x) 



H(x) 



^—i H'(a) x — a 



(1) 



where H(x) — Y\ ae ^{x — a) and H'(x) the derivative of 
H{x), has the property L r (a) = a\ r ' for all a G Q and is of 
the form L r (x) = x r ip(x n ) for some f(x) G F2[x]. 

Actually, Theorem A is an application of the following 
theorem which is also proved in PI . 

Theorem B: Let t be the error-correcting capacity and 
/ G F2 [xi, . . . , Xt] be a homogeneous symmetric function of 
degree r. Then there is a polynomial A e¥ 2 [x] such that 

f(z u ...,z t ) = A{{ Zl + ■■■ + zt) n ){zi + ■■■ + zt) r (2) 

for all ( Zl , . . . , zt) G T t = {(^ ,...,/?'", 0, . . . , 0) | < h < 
■ ■ ■ <l v < n — 1 and 1 < v < t}. 
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From Theorem A, the gap of the consecutive syndrome 
sequence can be filled in and then the binary cyclic codes 
generated by irreducible polynomials can be decoded by the 
BM algorithm. 

On the other hand, when the generator polynomial is not 
irreducible, the unknown syndromes cannot be expressed as 
a univariate function of any known syndrome. What we need 
here is a multivariate function. However, there is no, known to 
authors, multivariate interpolation formula (MVIF) over finite 
fields. In this paper, one such MVIF is developed based on 
the proposed Kronecker delta function over finite field. It is 
also proved that when applied to homogeneous symmetric 
rational functions on a certain set, the obtained polynomial 
has significant properties; namely, all the coefficients come 
from the base field and the degree of each term satisfies 
a congruence equation with modulus the code length. This 
is quite similar to that in the univariate case mentioned in 
Theorem B. Next, based on the proposed MVIF, the unknown 
syndrome representation method mentioned in H can be 
modified and then applied to any cyclic codes, which will 
be described in Section Hill 

Besides the unknown syndrome representation method, 
Orsini and Sala (8) presented an algebraic decoder based on 
the general error locator polynomial provided by the Grobner 
basis method, and for any cyclic code the correctable error 
patterns can be determined in one step. For those binary cyclic 
codes generated by irreducible polynomials, the general error 
locator polynomials can also be produced by applying the 
Lagrange interpolation formula |@). 

Furthermore, Orsini/Sala [9| and Augot et al. ifTUl , respec- 
tively, dealt with cyclic codes generated by reducible polyno- 
mials, in which the coefficients of error locator polynomials 
depending on more than one known syndromes. 

B. Motivation 

When decoding the binary quadratic residue code of length 
31 whose generator polynomial is a product of three irre- 
ducible factors, we obtained an explicit formula of Kronecker 
delta function and a multivariate interpolation formula over the 
binary field F2. The unknown syndrome representation for S3 
yielded from this MVIF is a binary polynomial and the degree 
of each term satisfies a congruent equation. Motivated by this 
example, a general result for arbitrary finite fields is presented 
and proved, and a one-step decoding algorithm is developed 
for any cyclic codes. 

C. Main Results 

In what follows, the original contributions of the paper are 
outlined as compared to the existing works. 

1) An explicit Kronecker delta function is proposed over 
finite field. For a G E = ¥ q m , if N = q m - 1, 

S a (x) — 1 — x N for a — 0, and S a (x) = — J2iLi( x / a Y> 
otherwise. 

2) A multivariate interpolation formula (MVIF) is devel- 
oped over finite field. Given M pairwise distinct vectors 
(x{,...,xf),...,(x 1 M ,...,xl I ) in E" = E x •■• x E 



and M arbitrary elements yi,..., um in E, the de- 
veloped MVIF yields a multivariate polynomial y = 
L(x x , . . . , x v ) over E so that for every i G {1, . . . , M} 
one has yi = L(x\,. . -,xV), 

3) We prove that, on a certain set, any homogeneous sym- 
metric rational function can be expressed as a polyno- 
mial function of some homogeneous symmetric rational 
functions. This polynomial can be provided from the 
proposed MVIF and has significant properties on both 
the coefficients and the degrees. 

4) For any cyclic code, the unknown syndromes can be 
expressed as polynomial functions of the known syn- 
dromes. The polynomials can also be provided by the 
presented MVIF and also have the same properties as 
mentioned in 3). 

5) By applying the unknown syndrome representation and 
the BM algorithm, a one-step decoding scheme is de- 
veloped for any cyclic codes. 



D. Organization 

In Section|nl all the theoretical results are presented, includ- 
ing the Kronecker delta function as well as the multivariate 
interpolation formula over finite fields and the relations among 
homogeneous symmetric rational functions. The decoder and 
the example based on the unknown syndrome representation 
are presented in Section [HI] The proofs of Theorems 2 and 4 
are provided in Appendix. 



II. Multivariate Interpolation Formula 

Throughout this paper, let n be a positive integer, q a prime 
power, m the multiplicative order of q modulo n. Also, let 
F = ¥ q and E = F 9 m be the base and extension fields of 
order q and q m , respectively. 



A. Kronecker Delta Function 

Since the set E* of nonzero elements in E forms a mul- 
tiplicative group, if N = |E*| = q m — 1, then 7^ = 1 for 
7 G E*, and if 7 G E*\{1} then 7 + 7 2 + • • • + 7^ = 



7 (1 



+ 7^-1) = 7 ((i - 7^)7(1 - 7)) = 0. Combining 



these two facts, one has an explicit Kronecker delta function 
on E. 

Definition 1. For a G E and N = |E| - 1, let S a (x) be the 
following function defined on E. 



S a (x) 



1 



„JV 



if a = 



It is easy to see that the value of S a (x) is 1 if x = a 
and otherwise. Therefore, 5 a (x) can be viewed as a finite 
field version of the Kronecker delta function on E. By using 
the function <5 Q (a;), a multivariate interpolation formula is 
developed over finite fields in the next subsection. 



3 



B. Multivariate Interpolation Formula over Finite Fields 

Based on the proposed Kronecker delta function in Defini- 
tion 1, one can develop the interpolation formula directly. 

Theorem 1. (multivariate interpolation formula over finite 
field) 

Let E be a finite field and N = E| — 1. For M a positive 
integer, if {(x{, . . . , x"), . . . , [x l M , . . . ,x v M )} is a set of M 
pairwise distinct w-tuples in E" = E x • • • x E and y\ , . . . , j/m 
are M arbitrary elements in E, then the multivariate polyno- 
mial given below 



y = L(x\...,x v ) 

has the property yi = L(x}, 
for a G E 



M v 

i=i j=i 

. , xf) for i = 1, . . . , M, where, 



Sa{x) 



1 — X 



N 



if a = 



The MVIF presented in Theorem 1 can be applied to yield 
polynomial relations among homogeneous symmetric rational 
functions and to represent the unknown syndromes in terms 
of known syndromes for arbitrary cyclic codes as mentioned 
in Theorems 2 and 4 below, respectively. 

C. Homogeneous Symmetric Rational Functions 

Since n divides q m — 1, there exist primitive nth roots of 
unity in E. Let f3 be one of them and let T t , t a positive integer, 
denote the following subset of E* defined as in 0, |4|: 



Tt = {(/?' 



,/3\0,...,0) | 0< h < 



< l v < n — 1, 
1 < v < t}. 



Recall that a rational function f(zi,...,z t ) over E is a fraction 
of two polynomials over F and denote by ¥[zi, . . . , z t ] and 
¥(zi, . . . , zt) the sets of all polynomials and rational func- 
tions, respectively, over F. A rational function f(z\, . . . , zt) 
is symmetric if /(^(l), ■ • • , ^r(t)) = f(zi,...,Zt) for any 
permutation tt of the set {l,...,t} and / is homogeneous 
of degree r if f{az\, . . . , azt) = a r f(z\, . . . , z t ) for any 
constant a G F. The fundamental theorem of symmetric 
polynomials says that any symmetric polynomials can be 
expressed as a polynomial function of the elementary sym- 
metric polynomials. For the homogeneous symmetric rational 
functions on the set T t , the following theorem indicates that 
one can express a homogeneous symmetric rational function 
f(z±, . . . , Zt) as a polynomial L(hi, ■ ■ ■ , h s ) over F in terms 
of some certain homogeneous symmetric rational functions 
hi(zi, . . . , Zt), . . . , h s (z%, . . . , z t ) and the degree of each term 
in L(hi, . . . , h s ) satisfies a certain condition. 

Theorem 2. (relation among homogeneous symmetric ra- 
tional functions) 

For positive integers s,ri,...,r s , let h\,...,h a G 
¥(zi,...,z t ) be homogeneous symmetric rational func- 
tions of degrees r%,...,r s , respectively. If all the vectors 
(hi(6), . . . ,h s (6)), for 9 G T t , are pairwise distinct, then 
for any homogeneous symmetric rational functions / G 



¥(zi,...,z t ) of degree r, there is a polynomial L G 
¥[x%, . . . ,x s ] of the form 



L(xi,. ..,x s ) 



E 



where J C {(it, ■ ■ ■ , i s ) | + • • • + r s i s = r (mod n)}, 
such that f(6) = L(hi(&), h s (6)) for any 6 G T t . 

Note that, the obtained multivariate polynomial 
L(xi,...,x s ) is over the base field F, that is, all the 
coefficients ,....j,'s are in F, and moreover, the exponents 
i%, . . . , i s of variables x/s in each term x 1 ^ ■ ■ ■ x\ s satisfy the 
congruence equation r\i\ + ■ ■ ■ + r s i s = r (mod n). It is easy 
to realize that Theorem B is a special case of Theorem 2 
with q = 2, s = 1, r\ = 1, and h(zi, . . . , z t ) = Z\ -\ + zt- 



D. Application to Decode Cyclic Codes 

Based on Theorems 1 and 2, algebraic decoders for arbitrary 
cyclic codes can be developed. As above, let n be a positive 
integer, q a prime power, m the multiplicative order of q 
modulo n, F = ¥ q , E = ¥ q m, and /? G E a primitive nth root 
of unity. Let C be a cyclic code of length n over F generated 
by g{x) — IliGSc ( x ~ wnere Sc is the defining set of C. 
Denote by rrii(x) the minimal polynomial of (3 l over ¥ q . If 
Rc is a subset of Sc such that g(x) = Y\ ieR nii(x) then call 
Rc a base set of Sc- Let t be the error-correcting capacity of 
C. An error pattern e(x) is called correctable if the number 
of nonzero terms in e(x) G E(a;), i.e., wt(e(x)) or the weight 
of e(x), is at most t. Let £ denote the set of all correctable 
error patterns. For i G {0, 1, . . . , n — 1}, call Si = e(f3 l ) the 
ith syndrome of the error pattern e(x). If i G Sc, Si can be 
calculated from the received word and is called the known 
syndrome of e(x); otherwise, for i ^ Sc, Si is called the 
unknown syndrome of e(x). 

Let s be the cardinality of Rc and let Rc = {r%, . . . , r s }; 
consider the set of s-tuples = {(e(/3 ri ),..., e(J3 rs )) \ 
e(x) G £ }. When C is the binary quadratic residue code 
of length 89, R c = {1,5,9,11} and it has been proved 
that there is a one-to-one correspondence between the sets £ 
and n = {(e(/3),e(/3 5 ),e(/3 9 ),e(/3 n )) | e(x) G £}. This is 
actually true for any cyclic code. 

Theorem 3. For any cyclic code C, there is a one-to-one 
correspondence between the set £ of correctable error patterns 
and the set £1 = {(e(^ ri ), ■ • ■ , e(/3 r =)) | e(x) G £}. 

Proof: It suffices to prove that for e\(x), e2(x) G £, one 
has e\(x) = e2(x) if and only if (ei (j3 ri ),..., e\(f5 rs )) = 
(e-i (P Tl ),■■■, e2(P Ts ))- Since the "necessary" part is trivial, 
we only prove the "sufficient" part. If (ei(/3*) : i G Rc) = 
(e 2 (j3 % ) : i G Rc), then ei(^) - e 2 (/3 l ) = for any i G 
Rc- Hence, for any i G Rc, mi{x) divides e\{x) — e 2 {x) 
which implies that g(x) divides e\(x)~~e2{x), or equivalently, 
e\(x) — e2(x) is a codeword, and then wt{e\{x) — e^ix)) > d. 
This leads to a contradiction because wt(e\(x) — 62(2;)) < 
wt(e 1 (x)) + wt(e 2 {x)) < 2t < d. ■ 

From now on, for convenience, T t will be replaced by 
T with t the error correcting capacity. Note that, in the 
binary case, i.e. q — 2, there is a one-to-one correspondence 
between the set £ of correctable error patterns and the set T, 
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with (/3' 1 , . . . , j8 " , 0, . . . , 0) G T corresponding to the error 
patterns e(x) = x' 1 + ■ ■ ■ + x' l! G £. 

Next, for i G {1, 2, . . . , n — 1} let Xj(zi, . . . , z t ) = z\ + 
• • • + zj! € F[zi, . . . , Zt] be the power symmetric polynomial 
of degree i. If 9 = (6 U . . . , 6 t ) G T, and e(x) G £ then 
Xi(9) = 9\ + ■ ■ ■ + 0\ gives the ith syndrome of e(x), 
i.e. Xi {6) = e(/3 l ) = S l . If R c = {n,...,r.}, then by 
Theorem 3 the power symmetric polynomials x n , . . . , x rs 
satisfy the condition of Theorem 2, i.e., (x ri (9), . . . , x Ts (9)) 
are all distinct for 9 G T. By Theorem 2 every homogeneous 
symmetric rational function can be expressed as a polynomial 
in terms of x ri , . . . , x rs . Since the unknown syndromes are 
also both homogeneous and symmetric, they can be expressed 
as polynomials in terms of the known syndromes x ri , . . . , x rs 
and this unknown syndrome representation is unified, i.e. 
independent of the occurred error pattern e(x). Hence, the 
decoding scheme in [4| based on the unified unknown syn- 
drome representation is valid for any binary cyclic codes. 
Furthermore, for arbitrary prime power q, the unified unknown 
syndrome representation method is also valid and will be 
shown in the following section. 

III. Applications to Unknown Syndrome 
Representation 

A. Unified Unknown Syndrome Representation 

If Rc = {ri, . . . , r s } is the base set of Sc mentioned above 
and S r is the unknown syndrome to be determined, then by 
Theorem 1 there is a polynomial L(x\, . . . ,x s ) over F such 
that S r = L(S ri , . . . , S rs ). When q = 2, since the syndromes 
can be viewed as homogeneous symmetric functions on the 
set T, by Theorem 2, the polynomial L(x\, . . . ,x s ) has the 
same form as mentioned in Theorem 2. For the general case, 
if q is a prime power, this fact is also true as shown below. 

Theorem 4. (unified unknown syndrome representation) 
The unknown syndrome S r can be expressed as a polynomial 
L(S ri , . . . , S rs ) in terms of S ri , . . . , S rs and the function L G 
F[xi, . . . , Xg] is of the following form 

L(x\, . . . , X s ) — ^ C-ii,...,i s X\ • • ■ X s 

(h,...,i s )£j 

where J C {(ii, ...,i s ) | r x ii H V r s i s = r s+1 (mod n)}. 

Since the proof of Theorem 4 is a slight modification of that 
of Theorem 2, only the sketch proof is given in Appendix. 

Now, with Theorem 4, the decoder based on the unknown 
syndrome representation in [4] can be modified and applied 
to any cyclic codes as shown in next subsection. 

B. Decoding Algorithm 

In this algorithm, we first calculate the known syndromes 
and use them to determine the needed unknown syndromes 
to obtain 2t consecutive syndrome sequence. Next, apply the 
BM algorithm to yield both the error locator polynomial a(x) 
and the error evaluator polynomial 51 (x), and then for v < t, 
the error pattern e(x) = e.^x^ 1 + ej 2 x^ 2 + ■ ■ ■ + ej^x^ can be 
determined by any appropriate root-position and -magnitude 
calculations. 



Input: r(x) 

1) Syndrome Calculation 

a) Known Syndrome Calculation 

Si =r(p i ),i G R c - 

b) Unknown Syndrome Calculation 

for all r < 2t, r <£ S c , 

S r = L r (S ri ,. . -,S r J, for {r 1: . . . ,r s } = R c , 

S r q = S^ . 

2) Error Locator and Evaluator Polynomials Calculation 

{a(x),Q(x)}=JFBMA(S 1 ,...,S 2 t). 

3) Root-Position and -Magnitude Calculation 
(Apply the Chien search and Forney algorithm) 

e{x) = 

for i from to n — 1 do 

if o-G0 _i ) = then e t = -Q(l3- l )/a'^/3" 1 ) 

Output: c(x) — r(x) — e(x) 



IFBMA is the abbreviation for inverse-free Berlekamp- 
Massey algorithm |5|. To illustrate this decoding algorithm, 
a workout example is given in next subsection. 

C. Decoding the (31, 16, 7) Binary Quadratic Residue Code 

Let E = F 2 5 be the finite field of order 32 with multiplica- 
tion modulus the primitive polynomial p(x) — 1 + x 2 + x 5 . 
Since (2 5 — 1)/31 = 1, if a is a root of p{x) then a 
is also a primitive 31st root of unity. Let C be the triple- 
error-correcting (31, 16, 7) binary QR code generated by the 
polynomial g(x) = 1 + x 3 + x & + x 9 + x 13 + x 14 + x 15 . Then, 
S c = {1, 2, 4, 5, 7, 8, 9, 10, 14, 16, 18, 19, 20, 25, 28} and the 
set Rc = {1,5,7} is a base set of Sc- Therefore, Si, S5, 
SV are known syndromes, and S3 is the needed unknown 
syndrome for the application of the BM algorithm. From 
Theorem 2, the unknown syndrome S3 = L(Si, S5, SV), 
where L(x, y, z) = J2(i j k)e.Ji ' zk an ^ listed in Table 
I, contains 307 hexadecimal ordered triples. For instance, the 
first triple (0,2, IE) in Ji, means the decimal ordered triple 
(0,2,30) which indicates y 2 z 30 is a term of L(x,y,z). Since 
C is a binary code, the only possible error value is one and 
there is no error-magnitude calculation. 

Now, suppose that the codeword c(x) = G C is 
transmitted and a corrupted codeword r(x) = x 3 + x 7 + x 20 is 
received, which indicates a weight-3 error e(x) = x 3 +x 7 +x 20 
occurred. Apply the decoding algorithm of Section IIII-BI as 
shown below: 

1) Calculate 2t consecutive syndromes: 

a) Compute known syndromes: 

51 = r(a) = a 3 + a 7 + a 20 = a 4 , 

5 2 = S 2 = a 8 , 

5 4 = S|=a 16 , 

5 5 = r(a 5 ) = (a 5 ) 3 + (a 5 ) 7 + (a 5 ) 20 = a 16 , 
S 7 = r(a 7 ) = (a 7 ) 3 + (a 7 ) 7 + (a 7 ) 20 = 0. 

b) Compute unknown syndromes: 

5 3 - L{S X , S 5 , S 7 ) = E ( ij lfc)eJa SI x Si x S 7 fc 
= S° x S| x S 30 + • • • + S 31 x S 5 13 x S° 

= a 27 , 

5 6 = S 2 = a 54 - a 23 . 
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2) Apply IFBMA to obtain the error locator polynomial: 
The inputs are the six consecutive syndromes 
Si,S 2 ,...,S 6 . While k = 0, set E°(x) = 1, 
A°{x) = 1, = 0, and 7C) = 1. Then apply IFBMA, 
as shown in Table [TTJ to obtain the error locator 
polynomial E^(x) = a 12 + a 16 x + a 3 x 2 + a n x 3 in 
the 6th procedure. 

3) Use Chien search or other method to find the error 
locations, and then correct the obtained errors. 

The three roots, a 11 , a 24 , a 2S , of (x) can be figured 
out and then the error locations are — 11 = 20, —24 = 7, 
and —28 = 3 (mod 31). Finally, the recovered codeword 
is c(x) = r{x) — e(x) = (x 3 + x 7 + x 20 ) — (a; 3 + x 7 + 
x 20 ) = which is as desired. 

IV. Applications to General Error Locator 
Polynomial 

The notion of general error locator polynomials was in- 
troduced by Orsini and Sala in O, (9), which is based on 
the following two facts: First, each coefficient of the error 
locator polynomial can be expressed as a function of the 
known syndromes. Next, there are one-to-one correspondences 
among the set of correctable error patterns, the set of known 
syndromes, and the set of error locator polynomials. The 
general error locator polynomial can then be defined in terms 
of the known syndromes and can be used to obtain the error 
locator polynomial without applying the Berlekamp-Massey 
method or any others. Orsini and Sala use the Grobner basis 
method to determine the general error locator polynomial for 
arbitrary cyclic codes by calculating each of its coefficients 
separately. 

In this paper, we apply the MVIF instead of the Grobner 
basis method to determine the general error locator polynomial 
also for arbitrary cyclic codes. Moreover, we prove that the 
degrees of each terms in the obtained polynomials satisfy a 
congruence condition. 

A. General Error Locator Polynomial 

For (9= (0i,..., t ) GT, let a(8, z) = (z-0 x ) ■ ■ -(z-d t ). 
Then a(6,z) = z* + dt-i^z*- 1 + ■■■ + ai(0)z + a o (0), 

where ^(0) = (-1)*~ 4 E fe , k± . . . fe ., Note that each 

coefficient cr, G W q [z%, . . . , Zt] is an elementary symmetric 
polynomial of homogeneous degree i, and by Theorem 2 can 
be expressed as a polynomial in terms of the syndrome polyno- 
mials xi(0),...,x s (0), i.e., di(0) = ai(xx(0),...,x s (8)) = 
di(x\, . . . ,x s ). The polynomial a(X, z) — a(x\, . . . , x s , z) 
obtained is called the general error locator polynomial of the 
cyclic code C. 

For any cyclic code, the coefficients Ui(x\, . . . ,x s ) defined 
above has the following form, as a consequence of Theorem 
2. 

Theorem 4. (Coefficients of general error locator polyno- 
mial) 

For 1 < i < t, the coefficient &i(xi, . . . ,x s ) of the general er- 
ror locator polynomial a(x\, . . . , x s , z) is in F g [xi, . . . , x s , z] 
and is of the form 

&j (x± , . . . j X s ) ^ Cii : ...,i s Xi ' ' X g 

(ii,...,i 3 )GJ 



where J is a subset of {(ii, . . . ,i s ) \ T\ix + ■ ■ ■ + r s i s = 
r s+1 (mod n)}. 

B. Decoding Algorithm 

To decode a corrupted data, we first calculate the known 
syndromes and use them to determine the coefficients of the 
general error locator polynomial. The obtained polynomial is 
actually the error locator polynomial. With some appropriate 
root position and magnitude calculations, the error can be 
pointed out and the original data can be recovered. 



Input: r(x) 

1) Syndrome Calculation 

for all i < 2t, S, = r(/3 l ) 

2) General Error Locator Polynomial Calculation 

a(x) = a(X,x), X = (S^, • • • , SfcJ, k\,...,k v G 

s c r\R. 

3) Root-Finding and Magnitude Calculation 

fl(x) ee S(x) x a(x) (mod x 2t ) 
e(x) = 

for i from to n — 1 do 

if <j(j3-' 1 ) = then e 4 = -VL((3- 1 ) / a' (fi^) 

Output: c(x) — r(x) — e(x) 



C. Decoding the (15,11,5) RS Code 

The decoding scheme based on the general error locator 
polynomial is illustrated in the following example. Let E = 
F 2 4 be the finite field of order 16 with the multiplication 
modulus the primitive polynomial p(x) = 1 + x + x 4 . If 
a is the root of p(x) then a is also a primitive 15th root 
of unity. Let C be the double-error-correcting (15, 11, 5) 
RS code generated by the polynomial g(x) — IliLi^ + 
a 1 ) = a 13 + a 6 x + a 3 x 2 + a 10 x 3 + x 4 . For the known 
syndromes S\, S2, S3, and S4, let the general error locator 
polynomial be a(X,x) = 1 + a 1 (X)x + a2(X)x 2 , which 
X = (S 1 ,S 2 ,S 3 ,S 4 ) and a x (X) = E^,*,-*)^ A.4A*X> 
a 2 (X) = T,(ij^ u )&A 3 x\xix^xf, where A 2 and A 3 are 
shown in Table [III] and Table [IV] 

Suppose a zero codeword c(x) = G C is transmitted, 
and the received corrupted codeword r(x) = a 6 x 2 + a 5 x 14 , 
which indicates that the error locations are 2 and 14, and the 
corresponding error values are a 5 and a 6 . Apply the decoding 
algorithm of Section IIV-BI as the following steps. 

1) Syndromes calculation: 

S t = r(a l ) = a 6 x (a l ) 2 +a 5 x (a 1 ) 14 where i=l to 4. 
Si = a 5 ,S 2 = a 12 ,S 3 = a 7 , and 5 4 = a 7 . 
S(x) = a 5 + a 12 x + a 7 x 2 + a 7 x 3 

2) General error locator polynomial calculation: 

^W = E (iJ , fe ,» )6 A 4 sixsixslxs^ 

x S ^ ~\~ ' ' ' ~\~ x tS'2 x 5g x 

— a 13 and 

^(x) = E (Wlfc tt ) g ^ ^ i5 x ^ x 5f x 

= S j x S2 x S x S ^ ~\~ ' ' ' ~\~ S-^ XiS'2 ^ ^4 
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= a. 

Therefore, a(X,x) = 1 + a x {X)x + a 2 (X)x 2 = 1 + 
a 13 x + ax 2 . 
3) Root-finding and magnitude calculation: 

The two polynomials <j(x) and Q.(x) satisfy the degree 
relation, deg Q(x) < deg cr(x) — 2, and £l(x) = 
S(x) x a(x) (mod x 4 ). Therefore, set Q(x) = a 6 + a 8 x 
and a'(x) = a 13 . One can find two roots of a(x) 
are a 1 and a 13 , which indicates that the locations of 
errors are 2 and 14. So e 2 = — Q(a 13 )/a 13 = a 6 
and ei4 = — ^(a 1 )/^ 13 = a 5 . The error polynomial is 
e(x) — a 6 x 2 + a 5 x u . Finally, the recovered codeword 
equals c(x) = r(x) —e(x) = 0, and finish the decoding 
algorithm. 

V. Concluding Remark 

The polynomials obtained from the proposed MVIF may 
contain quite large number of terms. However, in the decoding 
application, the situation may be mitigated. Actually, efforts 
have been made to reduce the number of terms in the unknown 
syndrome representation. One example is, by adding some 
extra points during the applications of MVIF to determine 
the unified unknown syndrome representation of the binary 
quadratic residue code of length 41, one may deduce the 
number of terms from 1355 in the original obtained poly- 
nomial to 1295 in a new obtained polynomial, and these two 
polynomials are equivalent, i.e. have the same values on the 
set £1 of known syndromes. This is just a first stage trial; that 
is, we just add the fewest possible points needed to yield a new 
polynomial in each MVIF calculation. We believe that some 
better polynomials, from the viewpoint of hardware design, 
will be obtained by adding more extra points in each MVIF 
calculation. 

Appendix 

A. Proof of Theorem 2 

Since all the vectors (h\(6), . . . , h s {6)), 9 G T, are pair- 
wise distinct, by Theorem 1, there is a polynomial L G 
E(xi,...,x s ) such that f(9) = . . . , h s (6)) for 

9 G T, and L(x\, . . . , x s ) is of the form 

s 

L( Xl ,...,x s ) = ^2f(^JlS hi{ ff ) (x i ), 

SeT i=1 

where S a (x) is the Kronecker delta function defined in Defi- 
nition 1. The polynomial L{x\, . . . , x s ) above can be written 

as 

L(xi, . . . , X s ) = ^ ] Cii,...,i s %i ' ' ' %'g i 

(»i,...,»«)e J 

where c iu ...^ s G E and J C {(ii,...,i s ) \ i\,...,i s 6 
{0, 1,2, ...}}. To prove the theorem, it suffices to show the 
following two claims: 

1) cn,...,i, S F for (ii, . . . ,i s ) G J, and 

2) J C {(ii, ...,i s ) | n«i H 1- r s i s = r(mod n)}. 



First, one defines some notations needed in the proof. For 

9 = (p u \...,p u -,0,...,0) 6 T and w £ {0, . . . , n — 1}, 
denote 

p w 6 = (p w \...,p w ",0,...,0), (3) 

where the set {w\, . . . , w v } equals {uT+w \ i = l,...,v}, 
Ui + w is the remainder of + w modulo n, and W\ < ■ ■ ■ < 
w v . Similarly, for w G {0, . . . , m — 1}, denote 

e qW = (^ i ,...,^«,o,...,o) ) (4) 

where {wi, . . . , w v } = {uiq w \ i = l,...,v} and w\ < • ■ • < 
w v . Note that since 7 9 ™ = 7 for 7 G E, if 1 ™ = if for 77 G T. 
Hence, there is a smallest positive integer d denoted by deg (9) 
such that 9 q = 9. Obviously, deg(9) is a divisor of m. Next, 
let "~" and "w" be two relations on T defined by 

~ $ if 0' = (3 W 9 (5) 

for some w G {0, . . . , n — 1} and 

9 » (9" if 9" = 9 qW (6) 

for some w G {0, . . . , m — 1}. Then, obviously, both "~" 
and "w" are equivalence relations on T. For 9 E T and d = 
deg(9), denote by (9) and [fl] the equivalence classes of 9 with 
respect to "~" and "w", respectively. Then one has 

(0) = I ^~ 0'} - {P w 9 \w = 0,...,n-l} (7) 

and 

[9] = 0" I (9 « 0"} = {^" |«; = 0,...,d- 1}. (8) 

Let i?i (resp. i? 2 ) be a complete set of representatives of "<~" 
(resp."w") in T, i.e. T can be expressed as a disjoint union 
of (9) (resp. [9]) for 9 G R\ (resp. ^ G R 2 ), or equivalently, 

7 1 = UeenM) ^ res P- T = UW^' Hence ' the polynomial 
X(a;i, . . . , x s ) can be expressed as either 

s 

6eRiae{8) i=1 

or 

s 

We will prove that for every 9 G R\ Claim (1) is true for the 
polynomial 

s 

L { ff ) (x U ...,X s )= ^ f( S )Yl d hi(a)(Xi) 

and for 9 G R2 Claim (2) is true for the polynomial 

s 

L [S] {X U ...,X S ) = ^ f(^)Yl S hi(a)(Xi), 

ae{3\ 1=1 

which will imply that both Claims (1) and (2) hold for the 
polynomial L(xi, . . . , x s ). 

For 9 = (^ Ul ,...,/3 U ",0,...,0) G T and a = 
F»§ = (/3^,...,/3-",0,...,0) G (9), if h(z u ...,zt) G 
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E 

w=0 



¥(zi,...,z t ) is homogeneous symmetric of degree r, 
then h(a) = h(/3 w 9) = h(/3 Wl , . . . , (3 W - , 0, . . . , 0) = 
/3 wr h(/3 Ul /3 U - , 0, 0) = (/3 w ) r h{9) = /3 rw h(9), since 

for each i, u>i = w + Uj (mod n) for some j. Similarly, for where 7 = f(0)((hi(9)) kl ■ ■ • (h s (ff)) k 'J- 1 G E. 

ff=#" G if fc(zi,. is a symmetric rational function Since 7^ = /(# J((/ii(0 9<i )) fcl ■ ■ ■ (h s (d qd )) k ->)- 1 = 

over F, then h(a) = h{9 qW ) = (h(9))^ . Hence, for either f(9)((h 1 (9)) kl ■ ■ ■ (h s (9)) k > )~ x = 7, this implies 7 belongs 

a G (9) or a G [9], hi(9) = if and only if hi(a) = for to the subfield of E of dimension d over F, and then 



i = l,...,s. 

Next, let Uq and U\ be subsets of {1, . . . , s] defined below: 



e;= 7 9 " g 



C/ = {i I - 0} and Ux = {1, . . . , s} \ C/ - 



Next, let A = (9). To prove Claim (2) 
is equivalent to show that all the coefficients 
E ae &f($)((hl(0)) kl ---(h s (6)) k *)- 1 in da are zero 
Let A be either (9) or [6]. The polynomial L A {x\, . . . ,x s ) except for those whose indices (kx, . . . , k s )'s satisfying the 
can be written as, from the definition of the Rronecker delta- congruent equation r\ki + • • • + r s k s = r (mod n). Now, 

since the rational functions f,hx, . . ■ ,h s are all homogeneous 
symmetric of degrees r, rx, . . . , r s , respectively, one has 



function, 
L A (xx, ...,x s 



E^ II S hi(S)( X i) II 6 hi{S){Xi) 
a£A ieU Q ieUi 

N / x k 



E/wlla-^IK-^E 



fc=l 

A? 



/it (a) 



E 



/(a) 



, ^(/»i(a))*--(^(a)) fc - 

ri-1 



E 



f{F°o) 



hi(a) 



Ild-^E/^IK-DE 

iet/ SeA ie(7i fc=i 



' o (fti0? u 'fl)) fc i---(/».(/3 to fl))*' 
P rw f{9) 

^ Q /3 r ^ (/ii(0)) fcl • • ■ /3 r > wk > {h s {9)) k s 



E 



71-1 



(9) 



f(0) 



-r s k s ) 



w=0 



where L^, . . . , x s ) = J2seA /(*) 11^ Ef=i(^ 

since n l6C /„( i - 



(a)> 



(hx(9))^---(h s ($)) k 



(hx(9)) 

^ ' pw(r—nki 



— -r.fc.J_ ( , ,, 



tu=0 



is a polynomial over the base Since /3 is a primitive nth root of unity, ^ A = unless 

if L* A (xx,-..,x s ) is a polynomial over F, \ is multiple of n. The formula in (HB equals zero, unless 



field 

then so does L A (xx, ■ ■ ■ ,x s ). Moreover, since n divides 
N = q m — 1, to show Claim (2) is true for L A (xi, . . . , x s ), 
it suffices to show the same assertion for L\(xi, . . . , x s ). 
Therefore, without loss of generality, we may assume 
L A (xx, ...,x s ) = L* A (xx, ...,x s ) and Ux = {1,.. .,s}. That 
is, 

s N , x k 

L A (xx,...,x s ) = J2f(Z)Ii(-VY, 



fci,...,fc. »=i v v ' 



M) s E E 



/(«) 



r — (rifci + • • • + r s fc s ) is a multiple of n, or equivalently, 
r\k\ + • • • + r s k s = r (mod ri). This completes the proof of 
Claim 2 and then the proof of Theorem 2. ■ 

B. Proof of Theorem 4 

The proof of Theorem 4 is almost the same as that of 
Theorem 2 with only some adjustments needed in the set- 
up step. First, replace the set T by the set of correctable 
syndromes fi = {e u j3 Ul + • • ■ + e Uv (3 Uv \ < U\ < ■ ■ ■ < 
u v < n — 1, e Ul ,...,e Uv G F, 1 < v < t}. Next, the 
relations "~" and "sa" on are simpler in this case; namely, 
for a, b G O, a ~ b if b = /3™a for w G {0, . . . , n - 1} and 



(hx(ct))^---(h s (a)) k » 



x kl ---x^. (10)a w 6 if 6= a?'" for w G {0, ...,d-l}, where d = de 5 (6,F; 

is the degree of b over F. Hence, if a = e Ul (3 Ul + ■ ■ - + e Uij /3 u " 

t n • m 1 ♦ A m a u * t and 6 = ^ a = e «x^ +Ul + ' ' ' + e u „/3 w+u Mhen 

To prove Claim (1), let A = [9\ and show that 

£ae[0] ^(") rL=i(Ma)r fcl G F for an Y index ( fc i> ■ • ■ > fc s )- 
Since deg(9) = d, [9] = {9 q \ w = 0, . . . , d — 1}, and since an( j 
f,hi,...,h s are symmetric rational functions, one has ^ ^ 

v /(g) ^ /(^) 



(/n(a)) fc i • • • (/i s (a)) fc = ^ {hx{9i w )) k ^ ■ ■ ■ {h s {9^)) k * 



e Ul P r{w+Ul) + ■ 
(3 rw a [r] . 



d-l 

E 



(f(o)r 



_Next, if b = of = e Ul /3 Ul 9™ + 



w =o((hx(e)) k ^ ■■■((h s (9)) k *)r 



= E 



b [r] 



e Ul /3 



ru\q 



e Uqj /3 u "i , then 



f(0) 



o\(M0)) k ^--(h s (9)y 



(aM) 9 ™. 



With these adjustments and some minor changes, the proof of 
Theorem 2 can be rewritten to obtain a proof of Theorem 4. 
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TABLE I 

List 307 triples (i, j, k) of Ji in hexadecimal representation 



0,2, IE 


0,3,10 


0,4,2 


0,5,13 


0,6,5 


0,8,8 


0,9,19 


0,A,B 


0,E,11 


0,10,14 


0,12,17 


0,13,9 


0,15,C 


0,1 A,4 


0,1B,15 


0,1E,A 


1,1,4 


1,8,1E 


1,9,10 


1,A,2 


1,C,5 


1,D,16 


1,E,8 


1,10,B 


1,14,11 


1,17,6 


1,18,17 


1,19,9 


1.1F.12 


2,0,9 


2,1, 1A 


2,2,C 


2,3,1D 


2,4,F 


2,5,1 


2,7,4 


2,A,18 


2,D,D 


2,F,10 


2,10,2 


2,11,13 


2,16,B 


2,18,E 


2,19,0 


2,1B,3 


2,1F,9 


3,0,1F 


3,1,11 


3,3,14 


3,6,9 


3,7,1A 


3,8,C 


3,D,4 


3,10,18 


3,16,2 


3,1A,8 


3,1C,B 


4,2,19 


4,5,E 


4,7,11 


4,A,6 


4,C,9 


4,10,F 


4,16,18 


4,18,1B 


4,19,D 


5,0,D 


5,1, IE 


5,2,10 


5,3,2 


5,5,5 


5,8,19 


5,10,6 


5,13,1A 


5,14,C 


5,17,1 


5,19,4 


5,1C,18 


5,1D,A 


6,0,4 


6,6,D 


6,8,10 


6,9,2 


6,A,13 


6,B,5 


6,C,16 


6,E,19 


6,113 


6,13,11 


6,16,6 


6,19,1A 


6,1 A,C 


6,1D,1 


6,1E,12 


7,0,1 A 


7,1,C 


7,4,1 


7,6,4 


7,8,7 


7,F,2 


7,13,8 


7,14,19 


7,18,0 


8,1,3 


8,2,14 


8,3,6 


8,4,17 


8,8,1D 


8,A,1 


8,D,15 


8,F,18 


8,14,10 


8,15,2 


8,17,5 


8,19,8 


8,1A,19 


8,1E,0 


8, IF, 11 


9,0,8 


9,1,19 


9,4,E 


9,8,14 


9,C,1A 


9,11,12 


9,12,4 


9,14,7 


9,15,18 


9, 16, A 


9,1C,13 


A,l,10 


A.3,13 


A,4,5 


A,7,19 


A,9,1C 


A,C,11 


A,D,3 


A, 12,1 A 


A,13,C 


A, 18,4 


A,1B,18 


B,0,15 


B,3,A 


B,4,1B 


B,7,10 


B,8,2 


B,10,E 


B,11,0 


B,18,1A 


B,1C,1 


C,1,1D 


C,D, 10 


C,13,19 


C,14,B 


C,1D,9 


D,9,l 


D,10,1B 


D,13,10 


E,4,0 


E,6,3 


E,8,6 


E,A,9 


E,C,C 


E,F,1 


E,10,12 


E,ll,4 


E,14,18 


E,19,10 


E,1E,8 


F,l,2 


F,A,0 


F,C,3 


10,0,7 


10,3,1B 


10,A,16 


10,B,8 


10,C,19 


10,D,B 


10,10,0 


10,13,14 


10,14,6 


11,3,12 


11,4,4 


11,5,15 


11,6,7 


11,8,A 


11,D,2 


11,E,13 


11,10,16 


11,13,B 


11,16,0 


11,17,11 


11,19,14 


11,1A,6 


12,2,17 


12,5,C 


12,C,7 


12,D,18 


12,E,A 


12,10,D 


12,12,10 


12,13,2 


12,14,13 


12,18,19 


12,19,B 


12,1 C,() 


12.1E.3 


13,0,B 


13.A.1A 


13,E,1 


13,13,18 


13,1D,8 


14,2,5 


14,4,8 


14,9,0 


14,C,14 


14,1 1,C 


14,16,4 


14,18,7 


14,19,18 


14,1E,10 


15,0,18 


15,2,1B 


15,8,5 


15,C,B 


15,F,0 


15,11,3 


15,15,9 


16,0,F 


16,1,1 


16,4,15 


16.8.1B 


16,C,2 


16,11,19 


16,15,0 


16,16,11 


16,18,14 


17,2,9 


17,7,1 


17,8,12 


17,9,4 


17,C,18 


17,11,10 


17,16,8 


18,1,E 


18,5,14 


18,6,6 


18,8,9 


18,9,1A 


18,D,1 


18,11,7 


18,13,A 


18,17,10 


18,18,2 


18,19,13 


19,0,13 


19,1,5 


19,3,8 


19,4,19 


19,5,B 


19,C,6 


19,10,C 


19,13,1 


19,14,12 


19,15,4 


19,19,A 


19, ID, 10 


19,1E,2 


1A,0,A 


1A,1,1B 


1A,6,13 


1A,8,16 


1A,E,0 


1 A, 10,3 


1A,12,6 


1A,19,1 


1A,1A,12 


IB, 0,1 


1B,1,12 


1B,2,4 


1B,5,18 


1B,6,A 


1B,A,10 


1B,B,2 


1C,0,17 


1C,2,1A 


1C,8,4 


1C,C,A 


1C,10,10 


1C,12,13 


IC.15,8 


1C,1A,0 


1C,1C,3 


1D,0,E 


1D,2,11 


1D,3,3 


1D,11,18 


1D,12,A 


1E,2,8 


1E,7,0 


1E,8,11 


1E,9,3 


1E,14,4 


1E,18,A 


1E,1C,10 


1F,4,2 


1F,8,8 


1F,D,() 













TABLE II 

Decoding Procedures While Applying IFBMA 



k 


AW 




A(»)(i) 




7 (fc) 





n.a. 


i 


1 





1 


1 


a 4 


l + a 4 x 


1 


1 


a 4 


2 





a 4 + a s x 


X 


1 


a 4 


3 


a 9 


a 8 + a 12 x + a 9 x 2 


4 i 8 
a + a°x 


2 


a 9 


4 





a^ + a 21 x + a ls x 2 


a 4 x + a s x 2 


2 


a 9 


5 


a 17 


a 26 + a 30 x + a 17 x 2 + a 25 x 3 


a 17 + a 21 x + a ls x 2 


3 


a 17 


6 







a 17 x + a 21 x 2 + a ls x 3 


3 


a 17 



TABLE III 
List 79 (i, j, k, u) elements of A2 



0,2,A,3 


0,4,6,5 


0,5,4,6 


0,6,2,7 


0,8,D,9 


0,9,B,A 


0,A,9,B 


0,B,7,C 


0,C,5,D 


0,D,3,E 


0,E,1,0 


0,E,1,F 


1,1,9,4 


1,5,1,8 


1,6,E,9 


1,D,0,1 


1,E,D,2 


2,3,2,8 


2,5,D,A 


2,C,E,2 


2,D,C,3 


2,E,A,4 


3,1,3,8 


3,4,C,B 


3,8,4,0 


3,9,2,1 


3,C,B,4 


4,3,B,C 


4,6,5,0 


4,B,A,5 


4,C,8,6 


4,E,4,8 


5,2,A,D 


5,5,4,1 


5,8,D,4 


5,A,9,6 


5,C,5,8 


6,1,9,E 


6,2,7,0 


6,8,A,6 


6,9,8,7 


6,A,6,8 


7,0,8,F 


7,1,6,1 


7,8,7,8 


8,6,8,8 


8,7,6,9 


8,8,4,A 


8,A,0,C 


8,E,7,1 


9,6,5,A 


9,8,1,C 


9,B,A,0 


9,C,8,1 


9,D,6,2 


A,2,A,8 


A,4,6,A 


A,5,4,B 


A,6,2,C 


A,C,5,3 


B,0,B,8 


B,4,3,C 


B,7,C,() 


B,8,A,1 


B,B,4,4 


C,3,2,D 


C,4,0,E 


C,A,3,5 


D,0,5,C 


D,2,1,E 


D,3,E,0 


D,4,C,1 


D,9,2,6 


E,0,2,E 


E,1,0,F 


E,8,l,7 


F,0,E,1 


F,7,(),8 


F,E,1,0 



























TABLE IV 
List 190 (i, j, k, u) elements of A3 



0,1,B,3 


0,2,9,4 


0,3,7,5 


0,6,1,8 


0,7,E,9 


0,D,2,0 


1,1,8,5 


1,3,4,7 


1,5,F,9 


1,6,D,A 


1,7,B,B 


1,8,9,C 


1,9,7,D 


1,A,5,E 


1,B,3,F 


1,F,A,4 


2,3,1,9 


2,6,A,C 


2,7,8,D 


2,8,6,E 


2,9,4,F 


2,A,2,1 


2,B,F,2 


2,F,7,6 


3,1,2,9 


3,2,F,A 


3,3,D,B 


3,4,B,C 


3,5,9,D 


3,6,7,E 


3,7,5,0 


3,7,5,F 


3,8,3,1 


3,A,E,3 


3,B,C,4 


3,C,A,5 


3,E,6,7 


3,F,4,8 


4,1,E,B 


4,4,83 


4,5,6,0 


4,5,6,F 


4,6,4,1 


4,7,2,2 


4,9,D,4 


4,B,9,6 


4,C,7,7 


4,F,1,A 


5,1,B,D 


5,2,9,E 


5,3,7,F 


5,4,5,1 


5,6,1,3 


5,7,E,4 


5,9,A,6 


5,C,4,9 


5,E,F,B 


5,F,D,C 


6,0,A,E 


6,1,8,0 


6,3,4,2 


6,5,F,4 


6,6,D,5 


6,7,B,6 


6,8,9,7 


6,A,5,9 


6,B,3,A 


6,D,E,C 


6,E,C,D 


6,F,A,E 


7,1,5,2 


7,2,3,3 


7,4,E,5 


7,5,C,6 


7,6,A,7 


7,8,6,9 


7,9,4,A 


7,A,2,B 


7,B,0,C 


7,B,F,C 


7,C,D,D 


7,D,B,E 


7,E,9,F 


7,F,7,1 


8,1,2,4 


8,3,D,6 


8,4,B,7 


8,7,5,A 


8,A,E,D 


8,B,C,E 


8,C,A,F 


8,F,4,3 


9,2,C,7 


9,3,A,8 


9,5,6,A 


9,7,2,C 


9,8,F,D 


9,9,D,E 


9,A,B,0 


9,A,B,F 


9,B,9,1 


9,D,5,3 


9,E,3,4 


9,F,1,5 


A,0,D,7 


A,3,7,A 


A,4,5,B 


A,7,E,E 


A,8,C,F 


A,A,8,2 


A,B,6,3 


A,D,2,5 


A,E,F,6 


A,F,D,7 


B,2,6,B 


B,4,2,D 


B,5,0,E 


B,5,F,E 


B,6,D,0 


B,6,D,F 


B,7,B,1 


B,8,9,2 


B,9,7,3 


B,B,3,5 


B,C,1,6 


B,D,E,7 


B,E,C,8 


B,F,A,9 


C,0,7,B 


C,1,5,C 


C,2,3,D 


C,3,1,E 


C,4,E,F 


C,6,A,2 


C,7,8,3 


C,8,6,4 


C,9,4,5 


C,B,F,7 


C,C,D,8 


C,D,B,9 


C,E,9,A 


C,F,7,B 


D,(),4,D 


D,2,0,F 


D,2,F,0 


D,2,F,F 


D,3,D,1 


D,4,B,2 


D,5,9,3 


D,6,7,4 


D,7,5,5 


D,9,l,7 


D,A,E,8 


D,B,C,9 


D,C,A,A 


D,D,8,B 


D,E,6,C 


D,F,4,D 


E,0,1,F 


E,2,C,2 


E,3,A,3 


E,4,8,4 


E,5,6,5 


E,6,4,6 


E,7,2,7 


E,8,0,8 


E,8,F,8 


E,9,D,9 


E,A,B,A 


E,B,9,B 


E,C,7,C 


E,D,5,D 


E,E,3,E 


E,F,1,F 


F,1,B,3 


F,2,9,4 


F,3,7,5 


F,4,5,6 


F,5,3,7 


F,7,E,9 


F,8,C,A 


F,9,A,B 


F,A,8,C 


F,B,6,D 


F,C,4,E 


F,D,2,() 


F,D,2,F 


F,E,0,1 


F,E,F,1 


F,F,D,2 













